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Abstract. We construct a no-arbitrage model of bond prices where the long 
bond is used as a numeraire. We develop bond prices and their dynamics 
without developing any model for the spot rate or forward rates. The model 
is arbitrage free and all nominal interest rates remain positive in the model. 
We give examples where our model does not have a spot rate; other examples 
include both spot and forward rates. 
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Models of term structure fall into two general categories: equilibrium models and 
no-arbitrage models. All the current models in each of these categories include spot 
rates and forward interest rates as part of the model structure (see [Mi] and [MR]). 
We introduce a no-arbitrage model for term structure that does not incorporate such 
rates. Although our model fits any reasonable data for the initial term structure 
exactly, some lattitude remains for a modeler to interpolate initial data so that 
interest rates may appear within the model. This option is a new feature for model 
calibration which might be exploited; for instance, implied interest rates may be 
developed from interest rate data to better fit the additional data. An example of 
this is given at the end of Section 6. 

We define a one-factor term structure model using a standard Brownian motion 
process which we denote by Bt- Our basic stochastic processes will be two well- 
known functions of Brownian motion: 

Mf will denote the simple exponential martingale 



Our bond prices will be defined in terms of the processes Mt and At, and the 
dynamics of these asset prices will be determined by the dynamics of these two 
processes. A scaling property of Brownian motion allows us to generalize our model 
in a straightforward way to include a volatility parameter (See Section 5). 

Our model places a special emphasis on a finite time horizon. Following the 
customary modeling approach, we calibrate our model at time t = with bond price 
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Mt = exp{Bt - -) 



and At will denote its time integral 
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data, the market price of discount bonds of various maturities. We use a standard 
notation (see [H JMj . for example) for these observed market prices, P{0,T), the 
parameter T being the maturity date. At maturity each bond is worth one monetary 
unit. The horizon date Th is the maturity date of the long bond. 

Definition 1.1. An horizon date Th is any fixed positive number for which there 
is an observable 'calibration value' for P(Q,Th). 

The market quote is called the long bond price at time zero and the model price 
for the long bond at any time t < T^ is denoted by P{t,Th). 

The stochastic process for the long bond is a crucial part of our model. This 
maturity sets the longest allowable maturity within the model as well as the time 
horizon for the evolution of the bond prices. Moreover, since the value at time t of 
our numeraire is P(t, Th), it is essential to define this process. 

However, one probabilistic feature of this model allows us to immediately de- 
scribe other maturity prices once the long bond process is set. This description 
involves a minimal amount of calibration data (whereas the long bond price in- 
volves a maximal amount of data) Therefore, we now describe, indirectly, the other 
maturity prices with the following 'out-of-order' definition: 

Definition 1.2. For any maturity date T < Th for which there is an observable 
'calibration value' for P(0,T), the stochastic process for the discounted T-bond price 
is given by 



Remark 1.3. The stochastic process defined in equation 11. 1\) is meaningful pre- 
cisely when the initial data satisfies P(0, T) > P{0,Th). Otherwise, a negative 
logarithm in the definition allows a collapse of the bond price to zero. This is a 
restriction on our model to reasonable calibration data. 

The initial discount bond prices must be monotone in the maturity date. This 
restriction, of course, does allow inversions of the yield curve, which is sometimes 
seen in real data. The restriction rules out only negative nominal forward rates 
within the calibration data. 

Equation (|l.ip will be used to define a T-maturity bond price once the long bond 
price process is chosen. Notice that the only data dependence within this stochastic 
process is the (constant) ratio of the initial bond prices. From Definition 11.21 one 
sees that this process has continuous, positive sample paths and the initial value of 
the process matches the ratio of the market bond prices at time 0. At this point 
we are able to illustrate our choice of the dynamics with the following proposition: 



(1.1) 



P{t,T) ^ /P(O^J 
P{t,Th) ■ \P{0,Th)) 



2 + At 1oe|P(0,T)/J=(0,T^)| 



Proposition 1.4. The process given by equation is a positive local martingale 
iffP{0,T)>P{Q,Th). 
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Proof. Let F ^ log{P{0,T)/P{0,Th)) and assume that F > 0. Equation can 
be written as 

P{t,T) _ / 2Mt \ 

We let 

y 2Af, 

and a basic calculation shows that 

(1.2) dY = YdB - ^Y^dt 

That is, 

f* 1 

Yt^F+ {YsdBs - -Ys^ds} 
Jo ^ 

This shows that is the exponent of a stochastic exponential (see [WH] or [KS] for 
a discussion of stochastic exponentials and their relation to the Girsanov theorem) . 
Hence, exp(yt) is a positive local martingale. In the case that F < 0, Yt has an 
initial negative value and explodes to — oo. In this case, the process in equation 
is not a positive local martingale. □ 

2. CALIBRATING THE LONG BOND 

Calibration data typically includes a finite sequence of bond prices correspond- 
ing to maturity dates Ti < T2 < ■■■ < Th. In fitting a model to this sequence, 
one frequently interpolates data and effectively determines an artificial bond price 
P(0,T) for all T < Th- For instance, one plots yields to maturity for the data at 
hand, and fits a yield curve through finitely many points. Then additional, artificial 
bond prices are determined from the completed yield curve. 

Some method of interpolating the observed bond prices must be carried out 
in order to define the stochastic process for the long bond price, P{t,Tfi)- Each 
method of interpolation gives arise to a different dynamics for the long bond. This 
in turn affects the dynamics of each shorter maturity (although it has no effect on 
the stochastic behavior of the discounted bonds in the model, as Proposition 11.41 
shows). 

Possible benefits of an interpolation scheme will be discussed in Section [6l For 
now, we assume that some method has produced a monotone decreasing function 
P(0, T) which passes through each observed bond price given by data. The simplest 
choice is the piece-wise linear curve through these data points. 

Definition 2.1. Given an horizon date Th > and a continuous, monotone de- 
creasing function T — > P{0,T) of 'calibration values' ,0 < T < Th such that 
P(0,0) = I, the stochastic process for the long bond price is given by 



4 



VICTOR GOODMAN AND KYOUNGHEE KIM 



(2.1) 



2 + At loglJ=(0,t)/J=(0,T^)| 



Remark 2.2. The stochastic process defined in equation \2.1]) has positive, con- 
tinuous sample paths on the time interval t < Th- Notice that the exponential base 
is a variable whose continuity properties depend explicitly on the continuity of ini- 
tial bond prices as a function of their maturity dates. Moreover, it is clear from 



Although one might choose the initial curve of maturities to be smooth, the pro- 
cess in equation h2.1\) does not inherit this smoothness because of the stochastic 
exponent. Nevertheless, as t ^ Th the bond price converges to one, as one would 
expect in a term structure model. 

Moreover, since the variable base is greater than one and the exponent is negative, 
the long bond price is always less than one before maturity. This is a feature that 
many models currently used in term structure modeling do not have (see [ MR, 
chapters 12 and 15] ). In our model, the yield to maturity is always positive. 

Definition \2.1\ is motivated by the discounted price expression in equation il) . 
// one formally sets the maturity parameter T equal to the time parameter t in hl.l\) 
and imposes the value P{t,t) — \, one obtains the reciprocal of the price given in 
Definition \2.1\ 

We use the formula appearing in Definition 11.21 and the long bond formula (|2.ip 
to define all the basic asset prices in the model: 

Definition 2.3. The stochastic process for each T-maturity bond price, T <Th is 



The time parameter t varies over the interval [0,r]. 



Proposition 2.4. For each maturity date T < Th, the T-maturity bond price has 
positive, continuous sample paths. Moreover, each price process exhibits continuous, 
limiting price behavior at its calibration time and at its maturity date. That is. 



equation 112. 1\) that 



\imP{t,Th) = P{Q,Th) 



(2.2) 




\miP{t,T) = P(0,T) 



and 



t~*T 



lim P{t, T) = 1 
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Proof. Arguments for the continuity behavior in the case T = Th were given in 
Remark 12.21 For T < Th, the product in equation p.2|) is clearly positive and 
continuous. First, we have 

hmP(,T).hmP(,T.)^.P(0,T) 



Secondly we observe from Definition 12.11 that 

is the expression for the reciprocal of the second factor in equation (|2.2p evaluated 
at f = T. This shows that 



lim Pit,T) = 1 



□ 



Remark 2.5. The proposition shows that each asset price evolves continuously 
from its initial value. In particular, if some Ti corresponds to an observed price at 
the time of calibration, then P{t, Ti) evolves from this price whatever interpolation 
scheme is employed to fill out the initial term structure (subject to the conditions 
of Definition \2. 

Of course, equation 112. 2(1 defines artificial bond processes for each interpolated 
value P{0,T). Remarkably, the time evolutions of these processes do not influence 
the processes P(t,Ti). This unusual feature of the model emphasizes the fact that 
instantaneous forward rates are not an essential feature of the model dynamics. 



In order to use techniques of stochastic analysis in developing our model, we 
establish a basic property of the dynamics given implicitly by equation (|2.2p . 



Definition 2.6. If a calibration scheme satisfies the conditions of Definition lKTi 
the initial price curve T P(0,T) defines the nominal forward rate curve F(t) 
where 

Fit)^iogiPio,t)/Pio,n)) 

Remark 2.7. F(t) is a continuous monotone decreasing function on [0, Th] and 
F{Th) ~ 0. Consequently, F(t) is a nonnegative function with finite variation. 



Proposition 2.8. 



an 



d 



(2.4) P{t,T) = exp( — -j^ — -r^ — 
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Proof. We may write equation ()2.1|) as 

\ogPit,Th) = -2F{t)Mt/{2 + AtF{t)) 

^ 2Mt 

2F{t)-^ + At 

and this establishes equation (|2.3[) . From equation (|2.2[) we have 

log Pit, T) = log P{t, Th) + 2F{T)Mt/{2 + AtF{T)) 
and equation (|2.4p is an immediate consequence. 



□ 



Theorem 2.9. iJac/i asset price process in the model defined by equations h2.1]] 
and 112. 2\) is a semimartingale. 

Each discounted asset price 

P{t,T) 

is a positive local martingale. 

If F{t) is a continuous, strictly decreasing nominal forward rate function with 
_F(0) = — log(P(0, Th)) then all nominal forward rates within the model are positive. 



Proof. Since the function F{t) has finite variation, the processes F(t)Mt and F{t)At 
are semimartingales (see [P, Chapter 2, Theorem 33]). It also follows from Theorem 
33 that any smooth function of these two semimartingales is also a semimartingale. 

Proposition 12.81 states both the long bond price and the T-maturity bond price 
are exponential functions of positive rational combinations of the processes. Hence, 
each asset price is a semimartingale. 

The local martingale property of the discounted assets was established in Propo- 
sition [TT2l 

If Ti < T2 the nominal forward interest rate aX t <Ti is 

log ( 



Tz-Ti °^P{t,T2)' 
Using Proposition 12.81 we obtain 



log P{t,Ti)- log P{t,T2 



2Mt 2Mf 



2F{Ti)-^+At 2F{T2)-^+At 



If F(Ti) > F(T2) this quantity is always strictly positive. That is, every nominal 
interest rate within the model is positive. 

Remark 2.10. The positivity of nominal rates also follows directly from equation 
il.l]) . One may use the equation twice and eliminate the extra term for log P{t,Th) 
to obtain the same formula for the nominal rate. 



□ 
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3. TRADING STRATEGIES AND ARBITRAGE 

We consider tame, self-financing trading strategies, similar to those introduced 
in [DH, Lemma 2]. 

Definition 3.1. // Ti < T2 • • • < Tn-i < Th are maturity dates, a trading strategy 
6 is a vector 

e= {61,62,..., On) 

where each function 6i is a bounded, progressively measurable function of t <Th 
and also satisfies 9i{t) — for t > Ti. 9i represents a position in the Ti-maturity 
bond; dn is the position in the long bond. Let 

Tg = inf{t<Th : e{t) = 0} 

A trading strategy must satisfy 6 = on the (possibly random) interval (T^i,T^]. 
A strategy 6 is said to be tame if the gains process 

G{t) := V / 6i{s)dP{s,Ti) 

is uniformly bounded below on the interval [0, Tg) when divided by the long bond 
price P{t,Th). That is, G{t A Tg)/ P{t AT0,Th) is bounded below by a non-random 
constant on the interval [0, Th] . 



Remark 3.2. Since a trading strategy is progressively measurable, the time Tg at 
which all positions are closed out is a stopping time. We restrict 6 so that 6 = 
for t > Tg because the model allows only bond trading. There is no money market 
in this model. 

Example. Let Ti < T2 < Th be maturity dates. Suppose that (Ti,i = 1,2, are 
stopping times with the property that CTj < Tj almost surely and cti < 02. We 
define another stopping time that controls the size of each discounted bond price, 
P{t,Ti)/P{t,Th). Fix C > and let 

(3.1) r = infO<T.: ^^^^^^-^^^ >C} 

The process in the definition of r equals log{P{t,Ti)/P{t,Th)) for t <Ti but this 
process is defined for all t <Th and it dominates all other logarithms of discounted 
prices with later maturity dates. Therefore, the other stopped process 

P{tAT,T2)/PitAT,Th) 

is also bounded. 

Let a and b be any constants and define a strategy where one holds a units of 
the Ti-maturity bond and b units of the T2-maturity bond until the random time 
ai. At this time one closes out the position in Ti -bonds and invests the proceeds 
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in the TJj-bond. The latter position is closed out at the random time CT2- Then 

ei{t) =al{4«,^} 

and 

{o-i<t<a-2} 

The gains process is 
(3.2) 

G{t)^ / adP{s,Ti)+ / bdP{s,T2) 
Jo Jo 

JaiAt P[(^l,-L2) 

= aP{ai A t, Ti) + bP{a2 A t, T2) + a pi'^'' [^(^2 A t, T2) - P{ai A t, T2)] 

However, if either a or 6 were negative, so that we have a short position in one 
bond, then the discounted gains process would not be bounded below. Neither 
of the discounted bond prices P{s,Ti)/ P{s,Th) is bounded so a constant short 
position, even up to a nonrandom time, is not a tame strategy. 

On the other hand, the strategy 

0{t)l{t<r} 

is tame. The indicator function modifies the original strategy so that when the 
time T is encountered, all positions are closed out. One can see from the formula 
for the gains process in equation (j3.2p that trading stopped at t effectively replaces 
each time ai by ai A t. All discounted quantities in the stopped version of equation 
(|3.2p are bounded. Therefore, the modified trading strategy is tame. 



Definition 3.3. A trading strategy 9 is self-financing iff the gains process G{t) in 
Definition \3.1\ satisfies 

9t ■ {Pit, Ti), . . . , Pit, Th)) = 00 • (P(0, Ti), . . . , P(0, Th)) + Git) 

fort<Tg. 

We will show that 

(3.3) 9t ■ iPit, Ti), . . . , Pit, n))/Pit, n) = 

9, ■ (P(0, Ti), . . . , P(0, r,))/P(0, Tu) + j2 f (^AS^J 

Equation (|3.3p is the well-known identity that relates self-financing portfolio values 
to stochastic integrals of the positions against the discounted bond prices. 

Consider the stopped version of equation p.3p where we stop trading at the time 
T defined in equation (|3.f p . The stochastic integral up to the first time f or r is 
effectively a stochastic integral up to time t with respect to the martingale processes 

Pis A r, T,) 
Pis A r, n) 

Each stopped process is a bounded local martingale according to Proposition If .41 
and hence is a martingale. Theorem 12. f of [Do] shows that a discounted portfolio 
value equals an integral with respect to the discounted asset prices under these 
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hypotheses. Therefore, the identity we seek holds on the event {t > t}. But t was 
defined with a parameter C. As we let C —^ oo, the event for which the identity 
holds expands to an event with probability one. This establishes equation (|3.3p for 
any self-financing trading strategy. 

Lemma 3.4. If is a tame, self-financing trading strategy, then the stopped, dis- 
counted portfolio value 

OtATe ■ {P{t A Te, Ti), . . . , P(t A Te, n))/P{t A Tg, T,,) 

is a supermartingale on the interval [0, Th] ■ 

Proof. Theorem 12 .91 states that each process P{t, Ti)/P{t, Th) is a local martingale. 
Since each position process 0i{t) is adapted and bounded, each stochastic integral 

is a local martingale. By assumption, the stopped, discounted gains process is 
bounded below and it is also a local martingale. It follows from Fatou's lemma 
that the discounted gains process, stopped at Tg, (and therefore the discounted 
portfolio value as well) is a supermartingale. 

□ 



Lemma 3.5. Let 

(3.4) Ilt^T^jPitAT,i,Th) 

denote a stopped, discounted portfolio value process, where all bond positions are 
closed out at some stopping time T^i < T^. Then the process also eguals the dis- 
counted value for bond trading which continues to the horizon date. 

If trading is extended by investing the monetary value T\ci in T^-maturity bonds 
at time Td and holding the bonds, then the process in \S.4^ is the discounted value 
throughout the interval [Q,Th]. 



Proof. At time a monetary amount IIj^^, may be exchanged for 

iiT^jp{Ta,n) 

units of the long bond. Thereafter, the portfolio value equals 

^''P{T,i,Th) 

Hence, the discounted portfolio value is also given by the process in equation p.4p 
for T,i<t< Th. 

□ 



Theorem 3.6. The collection of all tame, self-financing trading strategies for the 
asset prices given by Definitions \2.1\ and \2.S\ forms a system without arbitrage strate- 
gies. 
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Proof. Lemma 13.41 implies that each tame, seh- financing strategy produces a port- 
folio value so that its discounted version, when parametrized by i A Te is a super- 
martingale on the interval [0,r;i]. In addition. Lemma 13.51 implies that the process 
is a discounted portfolio value on this time interval. We apply [Do, Theorem 12.2] 
which states that if all discounted portfolio values are supermartingales there can 
be no (tame, extended) arbitrage strategy on the interval [0,r;i]. 

□ 

4. EXAMPLES OF INTEREST RATE DERIVATIVES 



Example 1. A Forward Contract on the Long Bond 

The obvious trading strategy to replicate a forward contract on a long bond 
allows us to illustrate a method for computing prices using local martingales. We 
denote the forward contract price for delivery at time T < Thhy f{t) . The exercise 
price required to receive a long bond is denoted by k. 

If the discounted gains process in equation (|3.3p were a martingale, then we 
would identify replicating portfolio values by computing the right hand expression 
in the 'identity' 

(4.1) = E[{P{T, n) - k)/P{T, n) I ] 

In our case the martingale assumption is incorrect, so we introduce stopping times. 
If one ignores the stopping time issue and computes a price using equation (|4.ip 
directly, one gets an incorrect answer (see appendix [B)) . 
Proposition 12.81 implies that 

r, = log[P(t,T)/P(,T,)].^^,^^,^ 

and for each n = 1, 2, . . . we define the random time r„ by 

(4.2) r„ = inf{i A T,, ; > n} 

It follows from an argument in the proof of Proposition 11.41 that 

exp(ytAr„) 

is a martingale. We then know that 

P{t/\Tn,T) 
P{tATn,Th) 

is a martingale. This suggests another contract that has a payoff similar to the 
forward, but may terminate early with a different terminal value. That is, we 
consider the conditional expectation 

(4.3) E[{P{T Kt,,,Th) - kP{T hT,,,T))/P{T KT,„n) \ Tt ] 

^P{TATn,n) ' ^ 

We think that the payoff is the long bond minus k shares of the T-bond and the 
payoff occurs at the smaller of the times T and r„ . This is useful because we can 
compute the conditional expected value in equation (j4.3p . It is 

P(iAT„,T) 



= 1 - At 



P{tAT^,Th) 
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Moreover, there is a tame trading strategy that produces the payoff value 

P(rAr„,r,0-«P(TAT„,T) 
It is 9{t) = {—K , 1), for t < r„ A T. Now the payoff value converges to 

with probability one as n — s- oo. And, the earlier price converges to 

Pit, Th) - KP{t, T) 

These limits are portfolio values because the trading strategy 'discovered' by com- 
puting in the stopped situation also converges io 6 — (— k ,1) ,t < T ] equation 
p.3p establishes the link between initial and terminal investment values via a trad- 
ing strategy, even in the limit as n — )■ cx3. 



Example 2. A Forward Contract on T'-maturity Bond 

In this example the delivery date is T and T < T' . The exercise price required 
to receive the bond is denoted by n. To pass to martingale computations we use 
the same sequence of stopping times as defined by equation (|4.2p in Example 1. 

As before, 

p(tAT„,r) 

exp(ytArJ - PI. rj.. 

P[t ATn,lh) 

is a martingale. We consider a contingent claim whose payoff is one share of the T'- 
bond minus k shares of the T-bond at the smaller of the times T and r„ . Equation 
()3.3p implies that a discounted portfolio, stopped at r„, is a martingale. Therefore, 
we consider 

S[(P(TAr„,r')-/tP(rAT„,T))/P(TAT„,Tft) | J^t ] 
(4.4) = E[{PiT ATn,T')/PiT ATn,Th) \ Tt\ 

- kE[P{T a Tn,T)/P{T A Tn,Th) I J't ] 

Each expected value is the conditional value of a martingale. The process with 
maturity date T is a martingale because the stopping time stops the local martingale 
so that it is a bounded process. Moreover, this process dominates the other local 
martingale point-wise because the process Yt is decreasing as a function of the 
maturity date. Therefore, the process with the T' maturity date is also a martingale. 
The conditional expectations in equation (14. 4|) equal 



P(tAT„,r) P(tAT„,T) 

K- 



P{thTn,Th) P{tATn,Th) 

so that an initial portfolio value for this contingent claim is 

P{tATn,T')~KP{tATn,T) 

We see that the trading strategy is 9{t) — (— k, 1,0) for t < T A t„. That is, some 
T-maturity bonds are sold short, a T'-bond is held, and no long bonds are involved. 
As n — > cxD the payoff, the initial portfolio value, and the trading strategy converge 
in a trivial sense. Equation (13. 3|) holds on the event {Yt < oo} which has probability 
one, therefore we see that the contract price equals 

P{t,T') ^ KP{t,T) 
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Example 3. A Caplet Price 

The payoff at time T' of a caplet for a LIBOR rate over the period [T, T'] can 
be written as [BGM, Section 3] 

Ct' - (T'-r)-i( ^ 



P{T,T') ' 

Here k = 1 + k{T' — T) and k is the payment 'cap'. Since the monetary payment 
(paid at the end of the LIBOR period) is set at the beginning of the period, it is 
simple to compute the market value, Ct, of the caplet at time T. The market value 
is determined by a fixed number of T'-maturity bonds. An investment consisting 
of Ct' many discount bonds at time T replicates the caplet payoff. So, we see that 

Ct = P(r,r')(r' - t)-\^^±^^ ~ ^)+ 

= (T' - T)-^{1 - kP{T, T'))+ 

As in examples 1 and 2, the stopping times r„ defined by equation (|4.2p allow 
each process 

P{tATn,f) 

to be a martingale ii T — T' or T. The first process is dominated by the process 
with maturity parameter T and both processes are bounded local martingales. 

Following the pattern used to compute prices in previous examples, we consider 
a different contingent claim that expires at the smaller of the times T and t„. The 
payoff , Cln, for this claim is 

CL = iP{T A T„, T) - nP{T A T„, T'))+ 

For convenience we omit the constant factor of (T' — T)^^ in defining this claim. 
Clearly, the limiting value of the payoff as n — > cx3 is the caplet value at time T 
when multiplied by the constant factor. 

We consider conditional expectations of the discounted payoff: 

(4.5) E[CI,JP[T hT,,,n)\Tt] 

The event D„ is defined by 

= {P[ThTn,T) > KP(rAr„,T')} 
The expected value in equation (|4.5p can be written as 

p(rAT„,r) p(TAT„,r) 

P[TATn,Th) P(TAT„,Th) 
Since each integrand contains a discounted bond factor, we may use Bayes' Rule 
[KS, Lemma 5.3] to compute each conditional expectation with a change of measure. 
The process 

P(^Ar„,f) ~ 
At = — -, exp — J I 
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is a positive martingale whose expected value equals one for either choice T = T' 
or T. At determines a probability measure Q„ where 

dQn X 

and Bayes' rule expresses the conditional expectation in equation (|4.5p as 
AtEgJ Id,, \^t]- i^KEq'J Id,, \ Tt ] 

^'''^ = ^P^M D. I ] - .^P^M A. I ] 

In Appendix we show that 

hm Q„ [ A. I ] - /C(T - t, log [^1^] , log k) 

2A'I 

where IC{s,x,y) = Pr{ _^ — ^— r— > ?/ or 2a::^^ < } 



(4.7) 



and 



(4.8) 



hm Q'^[D^\Tt]^)C'iT-t,log[^^],log^) 

n^ao P[t, 1 ) 

where /C'(s, a;, y) = Pr{-— > y} 



Therefore, we can take limits in equation (j4.6p to obtain the following price formula 
for the value of a caplet prior to the beginning of a LIBOR period: 



(4.9) 

Ct = r ip(t, r)/C(r, log [^^] , log n) S-'nPit, T')/C'(t, log [^^] , log n) 



where 5 = T' - T and t ^ T - t. 



5. THE VOLATILITY PARAMETER 



Each discounted bond price in Definition 11.21 is a positive local martingale. We 
will generalize this definition in a simple, but useful way. For any ct > it is obvious 
from Proposition II .41 that each process 

is a positive local martingale. That is, if we adjust the filtration by defining 

Gt = J^a^t 

then each process in equation (|5.ip is a local martingale. By speeding up (or 
slowing) the time evolution of asset prices with the factor ct we introduce a volatility 
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parameter within the model. This is analogous to the appearance of a volatility 
parameter in the widely used log normal model for a stock price. 

In fact, there is no need to rescale the time in order to include this desirable 
feature. The well known distributional equivalence of the processes {crS*} and 
{B„2^} allows us to modify Mt and At so that processes are effectively speeded up 
or slowed down. 



Definition 5.1. For a > let Af/'^' denote the simple exponential martingale 
and let A'f^ denote its time integral 

t 

M'f"'ds 

For a given horizon date Th as in Definition M.li the Th-maturity bond price of the 
general model is given by 

r -2M<''' 1 

2) P(t Th) ■— ( ^^^'^^ \ U+,t2a('') iog[p(o,t)/j>(o,r^)l-l 

^ ■ ^ V ' ■ Vp(o,r,0/ 

For a maturity date T < Th the T -maturity bond price of the general model is given 
by 

(5.3) P{t,T) ^ p(t,r/,)( ^^'-"'^^^ y^+-'-"i''-«'^<"-^'/^<"-^'^"^ 

for t < T. The model with cr = 1 will be referred to as the special model. 



Proposition 5.2. Each process 

P{t,T) 
P{t,Th) 

in the general model, determined by the asset prices given in Definition \5. 1[ has the 
same distribution as the process in equation i5.1]) . 



Proof. A simple time change using Brownian scaling shows that the exponent pro- 
cess 



2Mt 



2 + log(P(0,r)/P(0,T,)) 

of the general model discounted asset price has the same distribution (as a process) 
as the exponent 

2M,2t 

2 + A,2,log(P(0,T)/P(0,r,0) 
in equation (|5.1[) . □ 
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Corollary 5.3. The caplet price formula in the general model for a payment cap 
of size k on a LIBOR rate over the period [T, T'] is 

Ct =5-'P{t, T)/C(aV, log , log k) 

(5.4) ( 

- 6-'nP{t, T')/C'(aV, log [^|^] , log 

In this equation 5 — T' — T , 1 + 5k — n, and t = T — t. 



Proof. The proposition implies that the special model asset price process, parametrized 
by (T^t, has the same distributional properties as those of the general model. Then 
equation ()4.9|) . with the time rescaled, gives a replicating portfolio value for the 
general model. □ 



Corollary 5.4. An approximate value of the caplet price in the general model, valid 
for small values of a or small values of t = T — t, the time until the LIBOR rate 
is fixed, is given by 
(5.5) 

Ct ~ 5-^P{t, T)-^{bi) - 5-^KP{t, r')*(62) 

where 

_ \og[P{t,T)/P{t,r)]-\og. - /7wrp(, Typu T')^ 
a^\og[P{t,T)/P{t,T')] ± 2^^^^?>ynt,l)lt'[t,l )\ 

and 

$ denotes the standard normal distribution function. 



Proof. The function /C'(s, x, y) is defined in Example 3 of Section 4, and 

-2x-^+A„ 



IC'{a-r,x,y)^Pr{-^^^>y} 



which equals 



Pt{ — TT > y} 

The process Yr, where 

It — 



is determined by the SDE 

dYs = aYdB Y^ds 

2 

and its initial condition Yq = For small a we approximate Yt by Xt where X 
solves 

dX, = aYadB - —Y^ds 
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SO that 

2 

Xs = X + axBs —X s 

Then 

/C'((T^T, X, y) w Pr{x + ax^/rZ 2"'^^^ ~ ^"'^ 

where Z denotes a standard normal. Then 

/C' wPrja:: + ax^frZ 2"'^^^ ~ ^""^ 

=_Kr|(TV'^^ tt'^^ ^ j 

2 X 

=Pr{Z > + -V^x} 

<J^TX 2 

The other /C-function is approximately 



2A/*'^^ 



For this case we use the SDE 



dZ^ = (jxdB H x^ds 

2 

to approximate the process appearing in the definition of /C . 

^2 

/C w Pr{a; + ax^frZ + — x'^r > j/} 
We substitute the values 

Pit, T) 

X = log [ pi^^j.,^ ] and y = log[K] 

to obtain approximate values for the functions /C and /C' in the caplet formula in 
Corollary (|5.3p . This gives an approximate caplet price. □ 



6. FORWARD RATES AND THE SPOT RATE 



Definition 6.1. A term structure admits a family of a forward rate processes, 
r{t,T), if for each maturity date T < Th, 

(6.1) P{t,T) ^ exp{~ r{t,u)du) 

for t < T. 

Proposition 6.2. The special model admits forward rate processes if and only 
if the nominal forward rate curve F(t) satisfying Definition \ 2.6\ is an absolutely 
continuous function of t on the interval [0, Th] ■ 
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Moreover, if F{t) is absolutely continuous, so that 

i-Th 



F{t) = f{3)ds 

then for T < Th 

(6.2) r(t, T) = JJ^rj^y [ log[P(i, T)/P{t, T,,)]] ' 



Proof. Suppose that r{t, T) is a family of forward rates. If T < T/j we see from 
equation l|6.ip that 



P{t,T)/Pit,Th) = exp{ [ \{t,u)du) 
Jt 



Then it follows from Proposition 12 . 81 that 



r[t, u)du 



2F{T)-^+At 

This equation may be solved for the function F(T) and if the equation holds merely 
for some t on an event with positive probability, the equality implies that F{T) is 
the integral of an function. Therefore, it is necessary that the nominal forward 
rate curve be absolutely continuous. 

On the other hand, suppose that F{t) = J^'^ f{s)ds for some integrable function 
f{s). Then from Proposition 12 . 81 we obtain 

\og[p{t,T)/p{t,n)]^ 



2F(T)-i + At 
2Mt 



2[J^'^ f{s)ds]~^ + At 

which implies that 

log P{t, T) = ^ + log P{t, Th) 

2[J^-f{s)ds]-^+At 

This identity shows that the process log P{t, T) is an absolutely continuous function 
of its maturity parameter and that it is the integral of 

4Aft f{T) 



{2F{T)-^ + AtY F{TY 

That is, 

4M2 /(T) 



r{t,T) 



{2F{T)-^ + Atf MtF{Tf 

This establishes equation (|6.2p . The forward rates also represent the long bond. 
The formula above may be expressed as 

Af{T)Mt 



(6.3) r{t,T) = 



{2 + F{T)AtY 
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Then, noting that / ~ —F', we have 

r{t, u)du = 



At{2 + F{u)At) 
2Mt 



Th 



4Mt 



At At{2 + F{t)At) 
2F{t)Mt 



2 + F{t)At 

and this is the formula for — logP(t, Th) given in Proposition 



□ 



Corollary 6.3. // the nominal forward curve is absolutely continuous, then for 
t <T <Tii, each forward rate process r{t,T) satisfies the SDE 

(6.4) dr{t,T) = r{t,T)dB ~ r{t,T) r{t,u)du-dt 

Jt 

Proof. One may apply Ito's formula to the expression r(t, T) in equation (j6.3p and 
verify that the SDE holds. □ 

Corollary 6.4. Suppose that the nominal forward curve is continuously differen- 
tiable on [0, T/i). A continuous spot rate process r{t) may be defined as 

Proof. In |HJM| a spot rate r{t) is defined within an HJM model as lim^t r{t, T). 
If the function f{t) is continuous, it follows from equation 16.31 that 

limr(i, T) = lim '^^ ^ \ ! ^ n — -, f 

nt ^ ' ' Tit (2 + F(r)At)2 (2^F{t)Atf 

□ 

Remark 6.5. The existence of a spot rate process is merely an artifact of the 
special model. The model does not contain a money market asset; therefore the spot 
rate does not represent the rate of return on any investment. 

An assumption that the forward rate curve is smooth on the closed interval [0, T^] 
allows us to define the longest forward rate as the limit as T Th ofr(t,T). It 
follows from equation I16.3\) that 

r{t,n) ^ f{Th)Mt 

Notice that the longest forward rate process is a geometric Brownian motion. 



Example. The Cantor function, C(t), defined for < t < 1, is a well knonwn 
example of a non-absolutely continuous function [H]. C(t) is monotone increasing 
on the interval [0, 1], C(0) = 0, and C(l) = 1. We let Th — I and define a nominal 
forward curve function as 

Fit) = 1 - C(t) 

Then the special model determined by F(t) does not admit forward rate processes. 
In contrast, one choice for a smooth forward curve function is 

Fit) = ain - tf 
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where a and b are positive constants. The function f{t) = ab{Th — t)^^^ is integrable 
even in the case that 6 < 1 ; for this choice of 6 the initial instantaneous forward rate 
function, f{t), is convex. The special model determined by F(t) has the forward 
rates processes (see equation <!\6.3^ ] 

f. rj.. ^ 4abMt 

' {TH~TY-\2 + a{Th~TfAtf 
These rates blow up as the maturity parameter approaches the horizon date. This 
is an example with forward rate processes, but no forward rates to the horizon. 

Observations of instantaneous rate functions. The choice of a model horizon 
date may be a matter of convenience. On the other hand, one might set the model 
by selecting Th so that the corresponding long bond is actively traded and so that 
nearby maturities are also actively traded. This choice for Th would allow a large 
data set to be collected in a brief time period. One might use this to set the expected 
value of the longest forward rate 

It follows from the remark before the example that this expected value equals f{Th). 
Therefore, one might choose a nominal rate function F{t) which has a matching 
derivative at t = Th, fin) ^ ~F' (n). 

While observing the longest forward rate, one might also observe the spot rate. 
The stochastic portion of the spot rate in equation (|6.5p includes the same geometric 
Brownian motion that defines the stochastic behavior of the longest forward rate. 
From (|6.5p we obtain 

rU) ^^M^iog2[p(i Tft)l 



fin)f{t) 



log'[P{t,Th)] 



This implies that 
(6.6) 



rit,Th)F{ty 
r{t)r{t,Th) _f{Th)f{t) 



\oi[P{t,Th)] F{tY 
The identity in equation (|6.6p states that the deterministic function 

f{Th)f[t)/F{t)^ 

is determined by values of the long bond, the longest forward rate, and the spot 
rate. This is one example where a modeler might choose the function /(t), at least 
for t near zero, to approximate observations of r{t), r{t, Th), and P{t, Th). 



7. ORIGIN OF THE MODEL 

A risk-neutral measure in an HJM model [HJM, section 5] , [MR, equation 
(13.20)] produces the following dynamics for forward rate functions r{t,T): 

dr = a{t,T)J^ a{t,u)du ■ dt + a{t,T)dB 

The choice of 

a{t,T) ^ r{t,T) 
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is considered in [HJM , equation (45)]. The dynamics become 

dr — r ■ / r(i, u)du ■ dt + rdB 



Morton |Mo| found that this equation aUows interest rate explosions. Consequently, 
this model allows arbitrage. Explosions occur because of the quadratic behavior of 
the drift term. A positive initial value, combined with a volatility which prevents 
the rate from becoming negative, allows uncontrolled growth to infinity. 

Another dynamics for forward interest rates within an HJM model is given in 
[MR2, example 2.1]: 

dr = -(T{t,T) a{t,u)du-dt + a{t,T)dB 
Jt 

In this situation, when we specialize to cr(t, T) = r(t, T) we obtain 
dr = -r{t,T) r{t,u)du ■ dt + r{t,T)dB 



The appearance of the negative sign in the (quadratic) drift is due to a change of 
measure. This is the dynamics when a forward measure corresponding to the time 
Th is chosen. With such a measure, at least formally, each asset price discounted 
by the long bond price is a martingale. 

Solutions of this equation do not explode. And, conveniently, the solutions are 
explicit. If we state this equation in terms of the quantity 

Yt = r{t, u)du 



the equations reads 
This is the logarithm of the process we define in Definition 11.21 



dY ^ - W^dt + YdB 
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Appendix A. 

We identify two conditional probabilities that appear in Example 3, where the 
price formula for a caplet is obtained. In Example 3 we used stopping times r„ 
defined by equation (|4T2)) where <T <T' <Th. For i < T, the process 

~ P{tAT„,f) ~ 
A* = — -, exp — J I 
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is a positive martingale whose expected value equals one for either choice T = T' 
or T. At determines a probability measure Q„ where 



dP 

For each n — 1,2, . . . we let 



dQn T 

At 



_ P{TAr^,T) 



where k is a positive constant. 



Proposition A.l. 



lim Q„[ D„\Tt]^ Pr{ _^ * > logK or Ay-t > 2xt^} 

n^oo 2x^ — At-I 



where Xt 



P{t,T) 



(A.l) Pit,T') 
and 

lim Q'^[D^\Tt]^ Pr{ ^^^f-' > logn} 

n^oo 2Xt + AT-t 



Remark A. 2. The probability terms in the statement of the Proposition involve 
an independent process {Mg,As). That is, distribution functions of the random 
variables 

2M, 



2x-^±As 

are used to express the limit. 



Proof. We first consider the case T — T'. Then 

P(T AT„,Th) 

With respect to the measure Q'^, the process P{s, T)/P{s, T') is a Markov process. 
In order to see this, we use Proposition 12.81 to find that 

. P(g,r) , 2M, 2M, 

°^^P{s,T'y 2F(T)-i + A, 2F{T')-^+As 
which we write as 



Moreover, from Proposition 12 . 81 we also have 

= exp(y^^,^ - F(T')) 
Now both processes Yg and satisfy the SDE in equation (|1.2p . In particular, 

dYl^Y^dB^-^Yl^fds 
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This provides a very useful formula for A': 

A't = exp( Y'JBs - - {Y^fds) 

This shows that A' is a stochastic exponential determined by the bounded drift 
process y/ls<i-„. We apply the Girsanov theorem to see that A' adds the drift 
to the Brownian motion. That is, with respect to the measure Q'^ the process 

(A.3) B',=Bs- I " Y:,du 



JO 

is a standard Brownian motion. Now we derive an SDE for the process in equation 
(IA.2p . First, we have 

(A.4) d{Y - Y') = (r - Y')dB - _ {Y')^)ds 

since each process satisfies equation (|1.2p . However, under the process Bg is a 
semimartingale and we express the stochastic integral using equation (IA.3P : 

d{Y - Y') = (y - Y'){dB' + Y'ds) - - {Y'f)ds 
which holds path- wise for s < t„. We obtain 

d{Y - Y') ={Y - Y')dB' - -{-2YY' + 2(Y')'^ + Y^ - (Y'f)ds 
(A.5) 2 

={Y - Y')dB' - -{Y - Y'fds 

The autonomous form of this SDE shows that (Kj — Y^) is a Markov process. 
One immediate consequence is that 

Q'J D^\T^]^ Q'A Dn I log[||^^;Y^] ] 

However, a further consequence of equation (|A.5|) is that we can obtain an explicit 
description of Yg — Y^: 

Y -Y' ^ 

^ ' ^ " " 2{F{T) - F{T'))-^ + A'g 

In this equation, M' denotes the elementary exponential martingale for the Brow- 
nian motion B' and A' denotes its time integral. Now the event I?„ is 

= {>TAr„ - i"TAr„ > log 4 

and so we have 

a. \Tt]^Q'A YTAr„ - > l0g« I rtA.„ - i^Mr„ ] 

Since t„ is a stopping time for the process in equation (jA.6[) . we use the transi- 
tion probabilities for the process Y — Y' to compute the conditional probability. 
Denoting 

r 2M' 

by p{s,x,E), we apply [KS, Proposition 6.6, Chapter 2]: 

(A.7) Q'„[ Dn I .Ff ] = p(T A T„ - i A r„, Yt^r,, - Y^^. , [log k, oo)] 
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Equation (jA.7[) expresses the conditional probability in terms of a fixed function 
whose only dependence on n is through the stopping time t„. 
Since t„ ^ cx) almost surely as n — > oo, we have 

lim [Dn\Tt]^p{T^ t, Yt - F/, [log oo)) 

n — *oo 

This establishes the second limit in the Proposition. 

We next compute Qn[ Dn \ ] using a similar method. We have 

At = exp(rTAr„ - F{T)) 
and since the process Y satisfies the SDE given in equation (|1.2p . 

At = exp( / Y^dBs - - / Y^^ds) 
Jo ^ Jo 

so that the measure Q„ adds the drift Ys1s<t„ to the Brownian motion. With 
respect to Q„ the process 

=B,- I Y,du 



is a standard Brownian motion. From equation (jA.4[) we obtain 

d{Y - Y') =(y - Y'){dB + Yds) - _ {Y'f)ds 

={Y - Y')dB + i(y - Y' fds 

Again, the SDE for the difference process is autonomous, but now there is an 
explosion. The SDE has a positive quadratic drift term. An explicit solution is 

(A.8) - r; = '^-^ ^ 

2{F{T) - F{T'))-^ - As 
We denote the transition probabilities for this Markov process, 

2a;"i - As 

by q{s, X, E U {oo}). Then Proposition 6.6 of [KS] applies and we find that 

Qn[ Dn\Tt]^ q{T ATn-tATn, Yt^r^ - F/^^^ , [log K, Oo] ) 

As before, the only dependence on n is through the stopping times and we obtain 
lim Q„[ D„ I .Ft ] = q{T -t,Yt- Y;, [logK, cx)]) 

n^oo 

□ 



Pr{ e E or As > 2x^H, 



Remark A. 3. The exploding property of the process Ys — Y^ in equation iHT 



reveals that changing the numeraire can result in a singular measure. The measure 
Qn is defined in terms of the positive martingale 

P(tAT„,r) P{0,Th) 



A(t) 



PitAT^,Th) P(0,T) 
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Suppose that the equation 

dQ 1 PjO^Th) 

dP ■ P{T,Th) ' P(0,T) 

did define a probability measure Q. Such a measure would correspond, formally, to 
a change of numeraire where discounted asset prices of the form 

P{t,T) 

would be local martingales. However, for maturity dates S exceeding T , the price 
P(t, S) would be zero with positive Q-measure. Equation iA . ^) would give an ex- 
plicit formula for 

^ '^^P{t,T) 

assuming that 

B, =B,~ I Y„,du 



were a standard Brownian motion. But then the event — Ys = ~oo would corre- 
spond to {F{T) — F{T'))As > 2. Since Y^ — Yg is finite with probability one in the 
model, the measure Q could not be absolutely continuous with respect to the measure 
P. 



Appendix B. 



A naive use of martingale theory to compute prices for contingent claims may 
lead to incorrect answers. Here is one example where a forward contract price for 
the long bond can not be computed in the customary manner. The exercise price 
required to receive a long bond is denoted by k and the expiration date is T < T^. 

One might believe that the discounted contract price is given by 

(B.l) E[{P{T, Th) - k)/P{T, Th) I Tt ] 

^P{T,Th) ' ^ 

Theorem [2]9] asserts that the process P{t, T) / P{t, Th) is a positive local martingale, 
and it is well known that such processes are also supermartingales. Therefore, 

E[^^^^\Tt]<P{t,T)/P{t,Tu) 

In this case, the inequality is strict. To see this, suppose that equality holds with 
probability one. Then we have 

E[^^^^]^E[P{t,T)/P{t,n)] 

Proposition 12.81 allows us to write this identity as 



BARE TERM STRUCTURE 



25 



However, in [ GK, Corollary 4.2] we show that the expected value on the left hand 
side is strictly less than the other expected value. By using the Markov property 
of the process 

y 2M, 

* 2F{T)-^+At 
one can prove the stronger result that 

E[P{T,T)/P{T,Th) \Tt]< P{t,T)/P{t,Th) 

with probability one. Therefore, the computed price from equation (jB.ip is larger 
than the correct value 

P{t,Th) - KP{t,T) 
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